Universal Weak Lensing Distortion of Cosmological Correlation Functions 
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Gravitational lensing affects observed cosmological correlation functions because observed images 
do not coincide with true source locations. We treat this universal effect in a general way here, 
deriving a single formula that can be used to determine to what extent this effect distorts any 
correlation function. We then apply the general formula to the correlation functions of galaxies, of 
the 21-cm radiation field, and of the CMB. 
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I. INTRODUCTION 

Gravitational lensing has emerged as a powerful cosmo- 
logical tool. The spectacular images provided by strong 
lensing P, H, Q together with the robust parameter de- 
termination anticipated with weak lensing [J, [f| @, 0, H, 
0, EH EH are just two of the exciting developments in 
the field. When astronomers look for lensing, they have 
demonstrated that they can find it. Lensing can also con- 
taminate observations of correlation functions because 
the images of objects are displaced from the true source 
positions. For example, it has long been known that lens- 
ing smoothes out the spectrum of the cosmic microwave 
background (CMB) [3, 0, El- The effect of lensing 
on galaxy-galaxy correlations has been re-examined by 
several groups recently [H El El El El HI H3 • And 
a number of groups are anticipating the importance of 
lensing on future observations of neutral hydrogen via 
long wavelength measurements of the redshifted 21 cm 
line [13, [H, [2J] . 

In this work we present a unified treatment of the 
lensing contamination of correlation functions. First, 
we show that weak lensing contaminates any cosmologi- 
cal correlation function. Whereas previous studies have 
quantified its effect on some specific correlation functions 
(for instance, CMB temperature and polarization [l3| 
and baryon oscillations through the galaxy correlation 
function [l7l]). we focus here on the generality of the re- 
sult. This general framework allows us to understand 
when and for what correlation functions weak lensing 
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represents a sizeable or a negligible effect. By working 
in real space, we also hope to provide physical intuition 
for some of the properties of the lensing effects, which is 
harder to obtain from a pure fc-space calculation. Second, 
we consider the dependence of the lensing effect on the 
orientation of the source separation with respect to the 
line of sight. We also include the lensing- induced time 
delay as the longitudinal companion effect to the trans- 
verse lensing deflections. Finally, we show how the gen- 
eral framework presented here can be applied to different 
correlation functions: we retrieve previous results for the 
CMB, supplement and correct previous galaxy-galaxy re- 
sults [17] . and apply the formalism to the anisotropy in- 
duced by lensing of 21-cm radiation. 

The paper is organized as follows. In Sec. [ill a gen- 
eral formalism to calculate the weak lensing contribution 
to cosmological correlations is introduced along with a 
more physical understanding of the different quantities 
involved. In Sec. IIIII the general formalism developed is 
applied to three different cases: correlation of galaxies, of 
the cosmic microwave background and of the 21-cm radi- 
ation. Finally, Sec. IIVI contains some discussion together 
with concluding remarks. The calculations are relegated 
to the appendices. 



II. GENERAL FORMALISM 

Weak lensing, i.e. the effect of small potential differ- 
ences in the intervening space on the path of light, con- 
sists to first order of two effects: a transverse deflection 
of the photons from a straight path, and a time delay 
(Shapiro delay) along the path. Together, these effects 
result in a three-dimensional displacement of the appar- 
ent source position. 

Consider two physical observables A(x a ) and B(xb) 
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FIG. 1: Gravitational lensing bends photon trajectories. The 
sources, physically located at x a and xi are observed at x a +Ca 
and Xb + Ct (in direction n a and hb). a denotes the angle with 
the line of sight of the observed separation vector, f b s . 



and denote the observed values of A and B as A(x' a ) and 
B(x' b ). In this section we keep the nature of these cos- 
mological observables completely unspecified. The form 
of the results depends neither on the actual physical ob- 
servables that are being correlated nor on the nature of 
the source of the signal that is used to measure them. 
Since the geodesies light travels on are perturbed by the 
intervening distribution of matter, the measured values 
of the observables refer to physical points that are dis- 
placed with respect to the observed ones. 

Consider the arrangement depicted in Fig. [1] Since 
gravitational lensing perturbs the photons' trajectories, 
the two sources - located at x a and x b ~ are observed at 
S' a = %a+Ca and x' b = Xb + Cb- The observed distance r b s 
between the sources will in general differ from the true 
distance f. Whereas isotropy demands that the unlensed 
correlation function depends only on the distance r, the 
lensed correlation function depends on the angle a that r 
makes with the line of sight. Note that any difference in 
the line-of-sight angle of r and f obs is a higher order cor- 
rection, and we hence consider them to be equal for the 
perturbative approach adopted in this paper. Whenever 
the correlation function of two cosmological observables 
is measured, such measurement is subject to a modifica- 
tion arising because of the displacement of the apparent 
source positions x' a , x' b : 

(A(x' a )B(x' b )) = (A(x' a - Ca)B(x' b 00). (1) 

By Fourier transforming Eq. ([T]) and introducing the 
power spectrum of the two observables 

(A(fc a ,77 a )B*(fc & , 77 6 )) = {2tt) 3 5 (k a ~k b )PAB(ka,ria, r)b), 



it is straightforward to show that 
d 3 k 



(2) 



(A{x a )B{x b )) = 



ik-(x a —Xf,) 



(2tt) 3 
x ( e <£ -&-&), 



PAB(k,r] a ,T]b) 



assuming that the observables A, B are uncorrelated with 
the lensing deflection field (. Here r/A and r\B denote the 
conformal times at which light was emitted at the sources 
to reach us today. 

Equation §S§ allows us to express the effect that weak 
lensing has on the unlensed correlation function. In other 
words, it quantifies the effect that the matter distribution 
- responsible for lensing - has on the observed distribu- 
tion of the physical observables. Note that the weak lens- 
ing modification depends only on the difference £ Q — Cb of 
the deflections in both observables, not on the absolute 
magnitude of each. Assuming a flat Universe, as we do 
throughout, these distortions are given by the following 
integrals along the line of sight: 



Co, || 



2 / dx{Xa-x)^ 3 ^ X0a,X 



-2 / d X <S> 
h 



(4) 
(5) 



(3) 



Here the subscripts _L and || indicate directions trans- 
verse and parallel to the line of sight, respectively. In 
the transverse case, there are two such direction, indexed 
here with j = 1,2. The subscript a refers to the position 
of the observable A, fully specified by the direction 9 a and 
the comoving distance Xa- In the integral over x, the po- 
tential is to be evaluated along the line of sight so its 
argument is x(x) — [x@a, x) valid in the small angle limit 
we consider here. There is one important difference be- 
tween the parallel distortion £|| and the transverse distor- 
tion the latter has an extra factor of (xa — xY^ ■ The 
potential typically varies on small scales, so this factor 
is of order x/ r ^ 1- The transverse distortion therefore 
produces the dominant change in correlation functions. 

By neglecting all non-Gaussianity present in the lens- 
ing field, it is possible to evaluate Eq. ([3]) to all orders, 
thus generalizing the results obtained by Challinor and 
Lewis [l3[ in the case of CMB lensing. However, it is 
not clear to what extent this improves the accuracy, since 
non-Gaussian terms, induced by gravity during structure 
formation, appear in the expansion at all orders above the 
second and may play a significant role especially at low 
redshift and small scales. This issue can be resolved by 
tracing light rays through N-body simulations of cosmo- 
logical volumes (e.g., [25|]). In the case of CMB lensing, 
there are indications that deviations from the predictions 
of the calculation of [l3| which neglect non-Gaussianities 
in the deflection field begin to appear at I > 2000 in 
the deflection field. A detailed analysis of such a gen- 
eral treatment and of its range of validity, requiring ray- 
tracing measurement of the probability density function 
(PDF) of the displacements £, will be the focus of a forth- 
coming work. We therefore briefly outline this general 
formalism in appendix [Bl and for the rest of this work we 
will follow the perturbative approach. 

In order to proceed, we expand the exponential in 
Eq. Q to second order in the deflection. Note that 
all odd terms in this expansion vanish due to symme- 
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try, so that the expansion is good to fourth order. This 
is sufficient to quantify the weak lensing impact on the 
correlation function for sufficiently large separations. Ex- 
panding Eq. ([3]) above we then have 



(AB) « (AB) + 



d 2 (AB) d(AB) 



i' 3 



dr 2 



dr 



V 13 r 2 



(6) 



where r 



x\y denotes the observed separation be- 



tween the sources and the 3x3 distortion correlation 
tensor is defined as 

i^ij _ (CqCq) "f" (CfcCfc ) _ 



f T-D/2 -VJ 
T + D/2 

-v b o s 



(7) 



The second definition here of the matrix elements holds 
if we choose x a to lie along the z axis and both posi- 
tions to lie in the x — z plane, so that x a and Xb have 
coordinates x a = (0,0, Xa) and x b = (rj_,0, Xa + r \\) re- 
spectively. Note that due to the azimuthal symmetry of 
the lensing displacements around the line of sight, we can 
always choose such coordinates without loss of general- 
ity. In appendix[X]we derive expressions for the elements 
of the distortion tensor using the Limber approximation. 
The correlation of the distortions transverse to the line of 
sight decompose into a piece which is similar along both 
x and y directions 



T(x,r\\,r_i) = Ti(x,r\\) + T 2 (x,r±) 



with 



T2{x,r± 



dx 
o Jo 



dkk 3 P i> (k,x') 



d X ' 



2tt 

dkk 3 P$(k, X ') 



(8) 



(9) 



(x-x'f 



x (I - J {kr ±X '/x)) 



(10) 



and one which differs in the two transverse directions: 

"X 

'\2 



D( X) r±) = 2 I d X '(x-x'f 
Jo 

dkk 3 P$(k, X ') 



Mkr ±X '/x)-(U) 



Note that both T 2 and D, which depend only on the 
transverse distance between the two background objects, 
vanish in the limit that r±_ = 0. Since we will see that 
these two terms dominate the distortion of the correla- 
tion function, changes to the correlation functions have a 
characteristic dependence on the projected separation r± , 
peaking when r±_ = r. The variance of the displacements 
along the line of sight due to time delay is: 



S(x,r±) 



d x ' 



dkkP{k,x') 



x [1 - Mkr ± x'/X)} ■ 



(12) 



Finally, the displacement along the line of sight is slightly 
correlated with the transverse distortion; the relevant 
combination is 



V a + V b = 2r 



dx' 



dkk 2 P{k,x') 



Ji(kr ± x'/x)- 



(13) 

It is useful to estimate the order of magnitude of the 
corrections. First consider Eqs. (@J[5]). RMS of the 
perpendicular distortion is of order ^rmsX 2 / r ) where x 
is a typical cosmological distance, of order a Gpc, and 
r is a typical distance over which the potential varies, 
typically much smaller. The RMS of the parallel distor- 
tion is seen to be smaller by a factor of (r/x)- These 
estimates translate well into the respective variances as 
encoded in the functions T 2 and D. Both T 2 and D are 
seen from Eqs. (TlUl ITTj) to be of order k 3 P<s,(k)kx 3 since 
the line of sight variable x' is order x- The variance of 
the gravitational potential is roughly A|(fc) <~ fc 3 P$(fc), 
and the k integral picks out values of k ~ r^ 1 . So 
both T 2 and D are of order A|(l/r)(x 2 /r) 2 x (r/x), 
where A|(fc) is equal to $|_ M g ~ 10~ 9 on large scales 
(k < 0.01/i/Mpc) and suppressed on smaller scales due 
to the transfer function. The extra factor of r/x is 
the standard Limber suppression due to cancellation 
along the line of sight (only modes with k z small con- 
tribute to the variance). On large scales, the relevant 
dimensionless quantities T 2 /r 2 , D/r 2 thus are of order 
$| MS X 3 /r 3 - 10~ 6 (100 Mpc/h/r) 3 . They increase to- 
wards small scales, however not as quickly as r -3 , since 
the variance of the potential is suppressed by the trans- 
fer function towards smaller scales. In other words, due 
to the suppressed power on small scales, photons from 
both directions are deflected more and more coherently, 
and r b s is very close to r which is reflected by reduced 
values of T 2 and D. 

The first term in the transverse deflection correlation, 
T\ as defined in Eq. ©, represents the difference between 
the lensing experienced by A and B due to their different 
distances from us, so this part of the transverse lensing 
is suppressed by a factor of r 2 /x 2 and can usually be 
neglected. Similar estimates show that S and V a + V b 
are both smaller than the transverse variances T 2 and D 
by a factor of (r/x) 2 as expected from a cursory exami- 
nation of Eqs. I|4l5p . So we expect, and numerical work 
confirms this, that T 2 and D dominate the corrections to 
correlation functions, and the corrections will be of order 
A|(l/r)( X /r) 3 . 

Therefore, corrections to cosmological correlations due 
to weak lensing are small and given by: 



(AB)(r,r x ,x) 
T 2 (x,r±) 



, , , , ,2 d(AB) 

(AB r,x +T 2 X ,r x — 

r dr 



D(x,r±) \ r 
r 



2 dr 



1 d(AB) 
r dr 



(14) 



Two general features of this equation are worth point- 
ing out: first, if the correlation function is close to a 
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FIG. 2: Contour plot of T 2 /r 2 [Eq. JTDJl] determining the 
relative transverse displacement as a function of the source 
separation r and source redshift z (assumed equal for both 
sources, so r± = r). 



pure power-law, both derivative terms will be of order 
(AB)/r 2 . Hence, the relative lensing-induced effect on 
the correlation function will be given by T2/r 2 and D/r 2 , 
which are shown in Fig. [2] and [3J Second, in case of 
an oscillating correlation function, the lensing contribu- 
tion can be amplified by large values of d(AB)/dr and 
d 2 (AB) I dr 2 . Further, lensing will tend to smooth out the 
oscillations: at a local minimum of (AB), the first deriva- 
tive vanishes, while the second derivative is positive, so 
that the observed correlation is increased by lensing. At 
a local maximum, the opposite holds. This feature is 
already well-known in ^-space in case of the CMB. 

Equation (fT4"|) holds for any cosmological correlation 
function. It applies equally well to point-like sources and 
to diffuse backgrounds, to galaxies and QSO, to CMB 
and to the 21-cm radiation (albeit in the case of discrete 
sources magnification bias effects might provide the dom- 
inant distortion of the correlation function) . Given a par- 
ticular matter distribution, it is sufficient to evaluate the 
functions Ti and D once to be able to calculate the effect 
that weak lensing has on the correlation function of any 
cosmological observable. 

The functions T2/V 2 and D/r 2 are shown for a concor- 
dance ACDM cosmology with f2 m — 0.3 and Q\ = 0.7 
(which will be assumed throughout) in Figs. [5] and [3] We 
used the non- linear matter power spectrum from halof it 
(26| in the calculation. As the positions A and B move 
out to higher redshift, the lensing effects become larger 
since the longer path lengths lead to larger RMS deflec- 
tions. At fixed redshift, both T2 and D become larger 
as r gets bigger since the displacements of the path from 
the two points to us cease to be correlated and hence ex- 
perience independent deflections. However, as discussed 




FIG. 3: Same as Fig. [2] but for D/r 2 , where D is defined in 
Eq. (II). 



above, the relative effect of lensing distortions will gencr- 
ically be of order T^/r 2 , D/r 2 . These quantities decrease 
as r increases, as illustrated in Figs. [5] and [3J Hence, in 
the absence of features in the unlensed correlation func- 
tion considered, the effect of lensing is likely to be most 
important on small scales. Since the perturbative expan- 
sion treatment we are using is valid under the condition 
that {T 2 /r 2 ,D/r 2 ) < 1, Figs. [2] and [3j illustrate that the 
approximation used here is always applicable: it will yield 
a good approximation to the size of the lensing effect. If 
the effect is only marginally observable, this approxima- 
tion should be sufficient. If however the desired precision 
on lensing effect itself is high, as, e.g. in the case of CMB 
lensing in order to improve cosmological parameter con- 
straints, higher-order corrections will have to be taken 
into account. 

The value of the correlator of the longitudinal displace- 
ments S/r 2 = Z 33 /r 2 is shown in Fig.0J The longitudinal 
displacement effect (time delay) is clearly much smaller 
than the perpendicular one, as found earlier in [27l ]. 



III. APPLICATIONS 

We are now in a position to apply the above results to 
the correlation function of three different kinds of sources: 
galaxies, the CMB, and the 21-cm radiation background. 
The CMB is an angular measurement so all photons come 
from the same comoving distance meaning that ru = 
and the dependence is solely on r = rj_ . For galaxies and 
21cm, the orientation of the radius vector with respect 
to the line of sight is a free parameter, so the smoothing 
depends on r± even when r is fixed. 
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FIG. 4: Contour plot of S/r 2 [Eq. fl2])] determining the rel- 
ative longitudinal displacement as a function of the source 
separation r and source redshift z (assumed equal for both 
sources) . 



FIG. 5: Relative deviation from the unlensed correlation func- 
tion of the galaxy correlation function including effects of 
lensing deflections, near the baryon acoustic scale. Lensing 
smoothes the BAO feature in the galaxy correlation function 
at the percent level. 



A. Galaxies 

The effect of weak lensing on the galaxy correlation 
function has been taken into account in previous work 
in the context of the analysis of the impact that lensing 
has on the determination of the sound horizon of baryon 
oscillations [l7j . The ratio of the lensing induced term 
to the unlensed correlation function is plotted in Fig. [5] 
Lensing smoothes the BAO bump at the percent level. 
Note that since the lensing effect multiplies derivatives 
of the unlensed correlation function, the relative effect of 
lensing will be independent of any galaxy bias, in contrast 
to the magnification bias, whose relative effect is oc 1/6 
or cx l/b 2 , depending on redshift [l6|. We point out that 
Eqs. (6-8) of are recovered using Eq. (fl"4"|) above once 
the dependence of the kernels on the angle 9 is correctly 
taken into account in [l7l |. 

Fig. [5] shows the angular dependence of the correlation 
function for galaxies at redshift 3. The characteristic 
angular dependence illustrated in Fig. [6] may make this 
effect detectable. Note that due to the smallness of the 
time-delay effect (S/r 2 , Fig. 2]), the lensing contribution 
essentially vanishes for a — » 0. 



B. Cosmic Microwave Background 

The smoothing effect of lensing on the CMB power 
spectrum was computed long ago [121] in multipole space. 
The formalism established here allows for a simple cal- 
culation of this same effect in angular space. 

The angular temperature correlation function of the 



CMB is given by: 



w T t(o) = 

i=i 



2Z + 1 

4-7T 



C 1A (l)Pi(cosO), (15) 



where Pi denote the Legendre polynomials. Applying 
the approach presented here, we can calculate the lensing 
effect on the angular correlation function by evaluating 
Eq. (fT4"]l . Then, (AB)(r) stands for the temperature cor- 
relation function £,tt(t) = WTTfr/x*), where x* is the 
distance to the last scattering surface. For this, we need 
the first and second derivatives of £tt(0 which are given 
by: 



— sin 9 <Iwtt 
X* d cos 9 

~ COS 9 dWTT 



V 
A * 



dcos( 



sin 2 6 d 2 w T T 
? dcos9 2 



X 



(16) 



The derivatives of wtt with respect to cos 9 can be car- 
ried out using the Legendre polynomial relation: 



I 



1 



{xPi{x) - P^ix)) 



(17) 



Fig. [7] shows the lensed and unlensed CMB angular 
correlation functions in the region of the baryon acous- 
tic feature. One can discern a slight smoothing effect of 
lensing in real space, as pointed out after Eq. (fl"4"|) . Fig. [5] 
shows the difference between the lensed and unlensed 
correlation functions calculated in this approach (thick 
line). The unlensed C TT (l) were obtained from CAMB 
[28l |. In the figure we also show the correlation function 
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FIG. 6: Difference between the lensed and unlensed matter 
correlation function at z = 3 near the BAO feature. The 
three lines shown are for different orientations of the ra- 
dius vector with respect to the line of sight: perpendicular 
(solid/red), at 45° (dotted/green), and parallel to the line of 
sight (dashed/blue). The lensing contribution shows a strong 
dependence on the orientation of the separation vector. 

obtained from the lensed C TT (l) given by CAMB using 
Eq. (JT5J) (thin line). Clearly, the calculations in the two 
different approaches agree for 9 > 0.1° corresponding to 
r > 15Mpc//i. Note that results from N-body simula- 
tions [25j show deviations from the results of the CAMB 
code for i > 2000, roughly corresponding to 9 < 0.03°; 
however, it is not straightforward to convert the scale 
where deviations appear in multipole space to a corre- 
sponding angular scale in real space. 



C. 21-cm Surveys 

In principle, redshifted 21 cm radiation encodes infor- 
mation about the 3D distribution of neutral hydrogen in 
the universe. This distribution is sensitive to both inho- 
mogeneities in the matter and in the free electron density. 
It is not our purpose here to compute the complicated 
correlation function that results. Rather, we note that 
the real space calculation presented here is particularly 
simple to apply to any model of reionization (see Ref. [22j 
for a careful discussion of the complications that arise in 
Fourier space). 

For the purposes of this paper, we use the 21cm predic- 
tions for the "dark ages" before reionization calculated in 
[29l | . We use the 21cm spin temperature correlation co- 
efficients C TT (£) at z = 50, in the same way as outlined 




8 [deg] 

FIG. 7: Lensed (solid) and unlensed (dotted) CMB angular 
correlation functions. A slight smoothing effect washing out 
the baryon acoustic feature can be seen. 

above in the case of the CMB. Fig. [9] shows the relative 
lensing effect on the 21cm angular correlation function. 
It is rather small compared to that of the CMB due to 
the overall smoothness of the 21cm angular correlation 
function in the dark ages. Note that at lower redshifts, 
z < 12, we expect the 21cm correlation function to show 
a stronger baryonic signature (3(|, and hence expect a 
higher lensing effect possibly of importance to cosmolog- 
ical parameter constraints [3l|, [H, HH, HU • 

IV. DISCUSSION 

Gravitational lensing affects cosmological correlation 
functions in two ways. First, since the geodesies which 
photons travel on are sensitive to the distribution of mat- 
ter between source and observer and since the latter is lo- 
cally inhomogeneous thanks to structure formation, weak 
lensing acts by displacing the sources' observed (as op- 
posed to true) positions. The correlation function that 
is measured from observed data therefore also includes 
the contribution from these lensing deflections. Second, 
weak lensing also acts on observations by focusing or de- 
focusing the geodesies' congruences, thus altering the ob- 
served brightness of a source. This latter effect, known 
as magnification bias, is generally larger than the for- 
mer flq. I21I but it applies only to pointlike sources (see 
e.g. [l6l Il7l l20j for studies of the effect of magnifica- 
tion bias on two- and three-point correlation functions). 
Both these effects are small but potentially relevant in 
the present age of precision cosmology. In this work we 
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FIG. 8: Difference between the lensed and unlensed CMB 
angular correlation function wtt(6)- The thick line shows 
the calculation using the first-order approach presented here, 
while the thin line is the lensed wtt obtained from CAMB. 



FIG. 9: Relative lensing effect on the angular correlation func- 
tion of 21cm emission from the dark ages. We use the cross- 
correlation coefficients from [2S|. The redshift is z = 50, and 
the bandwidth chosen is 0.01 MHz. 



focused on the first of these effects and we derived a gen- 
eral perturbative formula that can be used to quantify 
it. 

The results of the present work depend only on the fol- 
lowing two assumptions: that photons travel on geodesies 
and that higher order corrections to the correlation func- 
tion (arising from higher order correlators of the lens- 
ing induced displacements) can safely be neglected. The 
first of this assumption has two consequences of oppo- 
site nature. The positive consequence is that the result 
obtained in this work is general and applies to any cos- 
mological correlation function, regardless of the nature 
of the source and of the physical observables that are 
being measured and correlated. Moreover, the correc- 
tion terms that appear in Eq. f| 14[) and that quantify the 
effect of lensing depend only on the power spectrum of 
matter: given a specific cosmological model they need 
to be evaluated only once and they can then be applied 
to any correlation function. As such, they also repre- 
sent a map that will tell whether weak lensing will have 
a relevant or an irrelevant role in the measurement of a 
given correlation function before the actual measurement 
is carried out. The negative consequence, on the other 
hand, is that this lensing-induced distortions represent 
a theoretical systematic that will always be present re- 
gardless of the precision of the instruments used to carry 
out the observation. In other words, even if "perfect ob- 
servations" free of any systematics could be carried out, 
this lensing-induced noise would still creep into the cor- 
relation function that would be calculated using those 



data. This lensing-induced modification of the correla- 
tion may be avoided only by reconstructing a map of 
the lensing potential (e.g., 35]). Another avenue to a 
possible disentanglement of this lensing distortion from 
the observed correlation function - at least for sources 
that are not confined to a fixed redshift - could be to 
exploit the angular dependence of the effect as seen in 
Fig. [S] Furthermore, since the lensing modification de- 
pends on the derivatives of the correlation function, cor- 
relation functions that are rapidly oscillating will in gen- 
eral be most affected by lensing, i.e. the features will be 
somewhat smoothed. Conversely, in quantities that are 
intrinsically uncorrelated, no correlation will be induced 
by lensing deflections. 

The second assumption entering the derivation of 
Eq. (fH| is that terms of third and higher order are dis- 
carded when the exponential of Eq. is expanded. As 
shown in Appendix [SI this assumption corresponds to 
Taylor expanding the coordinate dependence of the phys- 
ical observables and to retain only terms up to second 
order, which in turn corresponds to neglecting all con- 
tributions arising from the non-Gaussianity of the dis- 
placements' PDF. Despite the fact that under the same 
assumptions the general formalism outlined in Appendix 
|B]can be used to calculate the lensing distortion arising 
from the sum of all the terms appearing in the expo- 
nential of Eq. ([3]), it is however unclear to what extent 
this might really represent an improvement. If on one 
hand the sum of all the terms appearing in the exponen- 
tial may give important contributions on small scales, it 
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is also true that on such small scales departures from 
gaussianity of the displacements' PDF - induced by the 
non-linear growth of structure at low redshift - could 
play a relevant role. The range of applicability of such a 
non-perturbative method and the gain in precision that 
it would allow on small scales are the focus of a current 
project. 
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APPENDIX A: DERIVATION OF EQS. (|6 UT4)) 

1. Perturbative Approach 

Starting again from Eq. ^ it is possible to expand the 
exponential keeping only terms up to second order 

/AS]. 
_^_ e iH^-x b )p MkMb) 

X (l+ik- (Ca -(b) 

- ^[*- (&-&)][* •(£.-&)]), (Al) 

The zero-th order (in k) corresponds to the unlensed cor- 
relation function. Similarly, the first order term vanishes 
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because (() = 0. Let's then move to consider the second 
term, which we'll denote as (AB) 2 . We can rewrite it as 



(AB) 2 = - 



d 3 k 



M*.-*)p AB (k,Tla,ri b ) 



(CCD 

d 3 k 



^ik-(x a — x b ) 



PAB(k,rja,r]b)kikj (A2) 



where we have used the fact that the displacement cor- 
relators do not depend on the integration variable k to 
pull them out of the integral. The above integral can be 
rewritten as 

d 3 k - 

T e lkr P A B(k,r) a ,'q b )k t k : j 



(2tt)3 

d 3 k 1 d 2 
(2ir) 3 i 2 dridrj 
1 d 2 f d 3 k 



Ak-r 



PAB(k,7] a ,7] b ) 



i 2 drjdj 



Ak-r 



(2tt) 3 



PAB(k,r) a ,r)b) 



d 2 



(AB). 



dridrj 
We can then notice that 

d 2 {AB) _ d 2 (AB) dr dr d(AB) d 2 r 



(A3) 



dridrj dr 2 dri drj 

nrj d 2 (AB) 

l d(AB) 

+ 1 °- 

r dr V " r 



dr dridr.; 



(A4) 



Putting all the pieces together we then get 



(AB), 



{CaCD + iClQ) 



(CCD 



nrj d 2 (AB) 1 d(AB) 



dr 2 



dr 



(A5) 



which is exactly the result of Eq. ([6]) once we define the 
displacement correlator as 



(ca> + (act) 



(A6) 



Finally, let's notice that the same result can be ob- 
tained in a somewhat more intuitive way simply by Tay- 
lor expanding the coordinate dependence of the observ- 
ables as 

A(x a ) = A(x a + Ca) « + £A°, + \Q(IA% , (A7) 

where we use the shorthand notation A^ = dA/dx l | s _ s , 
and where the "a" index appearing on the displacement 
and on the observables A a specifies that these quanti- 
ties refer to the physical point x a . 



2. Decomposition of the Displacements' Correlator 

Consider a perturbed flat FRW spacetime with metric 

ds 2 = a 2 (ri) [(-1 - 2tf)dr? 2 + % (1 + 2$)dx i dx j ] . 

(A8) 

The lensing induced deflection of a source at distance 
Xa is given by the following integrals along the (unper- 
turbed) line of sight d, [36| 



rXo 

C ± = 2 / dxixo - x)Vi$ 
Jo 



rXo f d 3 k 

= 2i J d X (xo-x) J j^ki$e ik - s M, (A9) 
rXo 

Ca,\\ = -2 / d X <S> 
Jo 

= - »jf < Aio > 

With the help of the Limber approximation, it is straight- 
forward to show that the correlators are given by 

/■min{ x °,X?} 

(Co,||C&,||) = 4 / d X 
Jo 

dk 

I —kJo{k X e)P^(k,x), 



2ir 



(All) 



(CLlC 



'.'■11/ 



4« / d X {X° a ~ X) 

Jo 

I £^y ± e^^P^[k ±1 z(x)] (A12) 



(27T) 2 
r mm{x a ,X°} 

(£,±&±) = 4/ dx(xl-x)(x° b -x) 



d 2 k ± 
(2n) 2 



feifcl e^xp*^, z( x )], (A13) 



where in Eqs. (|A12[ [AT3j) i,j — 1, 2 denote the two com- 
ponents perpendicular to the line of sight. Looking at 
the structure of Eqs. (|A1HIA13|) above it is possible to 
notice that these quantities transform as a scalar, a vec- 
tor and a symmetric tensor with respect to rotations in a 
plane perpendicular to the line of sight. We can therefore 
define right away the scalar S a b as 

Sab = {Q,\\a,\\), (A14) 

and can simplify even more Eqs. (|AT2l rAT3|) by extract- 
ing the components of the vector and of the tensor. In 
the case of the vector, for instance, it is possible to notice 
that the vector component is aligned along the displace- 
ment vector (C^ j_Cb,||) ~ We can then define 

Va = ^(CLCm) (A15) 

r mm{ X ° a ,X°} 

dx(Xa - X) 

/k 2 dk 
—P 9 [k,z( X )]Ji(k X 0). (AI6) 
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Similarly it is possible to decompose the tensor part into 
its trace and its off diagonal traceless symmetric part (cfr. 
[13]). Letting 



(A17) 



where R? is the symmetric traceless tensor defined by 



(A18) 



and then using 



(A19) 



\T af ,-^) (A20) 



together with the Limber approximation, it is possible 
to g et to the following expressions for T a b and D a b (cfr. 



Tab = 4 / dx(Xa"X)(x2-x) 
Jo 

k 3 dk 



x 



-J (k X 6)P^k,x), 



4tt 

r-min{ X °, X ?} 

£U = 4 / dx(x°-x)(x°-x) 





k 3 dk 

2tt 



J 2 {k X 6)P*{k, X ), 



(A21) 



(A22) 



Taking the limit 6> -> of Eqs. (jATTl lATHl |A122|) it is 
possible to obtain the equivalent of the above expres- 
sions for the single source case. In particular, given that 
Ji(0) = J 2 (0) = 0, it is straightforward to notice that in 
this case V a = and D aa = 0. We can then forget about 
the labels for D a b and from here on simply identifying it 
with D. 

Let's notice that up to this point only the Limber ap- 
proximation has entered the above derivation. It is then 
possible to proceed further by defining 



T = — (T aa + Tbb — 2T a b) 



(A23) 



Defining the following shorthand notation for the sake of 
brevity, 

^k 3 P*(k, x ) 7 (A24) 

/rlk 
— k 3 P^k lX )Mkx6), (A25) 



we can then rewrite Eq. (|A23[) above as 



T = 2 



Xa j-Xb 

dx(Xa ~ x) 2 Vx) + / dxiXb - x) 2 Io(x) 
o Jo 



dxiXb ~ x)(Xa - x) Ij(x) 



dx[(Xa~x) 2 + (Xb-x) 2 ]lo(x) 



dxiXb ~ xfhix) 



- 2 / dx{Xa - x){Xb - x)Ij(x) 
la 



dx(Xb ~ Xa) 2 Io(x) 



Xb 



dxiXb ~ xfloix) 



2 / dx{Xa - x)(Xb - X)[I.J{X) ~ MX)} 
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(A26) 



where without loss of generality we assumed that Xa < 
Xb- Equation (|A26[) above is exact. We can proceed fur- 
ther by noting that Xb~Xa = f\\ and that the integration 
limits of the second integral extend over ru . We then have 



Xa rXb 

dxiXb - Xa) 2 Io(x) + / dxiXb - x) 2 Io{x) 



Xa fXb < _ \2 

d X Io(x)+ dxI (x V X r 



Xa '|| 

k a P*(*,x)=2Ti(x a ,r||), (A27) 



where to get to the last line we used the fact that 
while the first integral scales as r^Xa, the second one 
scales only as ru and it therefore provides a subdomi- 
nant contribution. Finally, the last term of Eq. (|A26D 
gives T2( Xa ,r±) provided that we make the approxima- 
tion (x - Xb) ~ (X - Xa)- Then 

Xa 

dx(Xa - x)(Xfe - X)[lj(x) ~ Io(x)] 

Xa r JI, 

dx(Xa~x) 2 / -fc 3 P*(fc,x) [l-Jo(k X 0)} 
= T 2 (Xa,r±). (A28) 
Putting all the pieces together we have then shown that 
T = T 1 ( X a,r ll ) + T 2 (xa,rx)- (A29) 

Finally to recover the matrix decomposition of the Z v 
correlator we can appeal to the simmetry of the lensing 
displacements £. Notice in fact that from Eqs. (|A91IA10"|) 
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the displacements are clearly invariant with respect to 
a rotation around the line of sight. We can then arbi- 
trarily pick the coordinate system so that the sources lie 
in the x — z plane, with the z axis directed along line 
of sight to x a . The coordinates of the sources are then 
x a = (0,0, x a ) and Xb = (rj_, 0, Xa + »"||) respectively, 
and since the displacement vector r makes an angle a 
with the line of sight we have that rj_=r sin(a) and 
rii = r cos(a). From here it is then straightforward to 
obtain the expression for as in the second line of 
Eq. 0. 



3. Derivation of Eq. ([14]) 

Finally notice that the correction term appearing in 
Eq. ^ due to lensing can be recast in the form 



(AB), 



Zn 



nr 



d 2 (AB) d(AB) 
dr 2 dr 



d (\d{AB)\ 
r 



dr \r dr J 



d{AB) Hj 

r — d J 

dr 



(A30) 



Then, using Eq. ([7]) for Z' 1 ^ , which is valid with this choice 
of coordinates, we get 



T 



D 



r1 (v a + v b ) r ^ 



1 d(AB) 
r dr 



(2T + S) , 



(A31) 



which reduces to Eq. (TH)) once we notice that T2 and D 
provide the dominant contributions. 



APPENDIX B: THE EXACT BUT MORE 
RESTRICTIVE TREATEMENT 

Following the approach of [l3|, it is possible to obtain 
an "exact" expression for the exponential appearing in 
Eq. ^ under the restrictive condition that the quantity 
ki{C a — Q) is gaussian distributed. If this condition holds 
then we can use the fact that 



(e ly ) 



-<y 2 >/2 



(Bl) 



to rewrite Eq. §5§ as 



(AB) = 



d 3 k 



(2tt)3 

X e -3([fe-(Ca-Ci.)] 2 > _ 



PAB(k,r] ai r) b ) 



(B2) 



We therefore need to calculate the value of the correlator 

- i 2\ 



k ■ (Ca - Cb) 



{(ICi) + {QCD - {CCD 



(cici) 



(B3) 



It is possible to proceed further by taking into account 
the fact that the lensing-induced displacements are char- 
acterized by different expressions that depend on the di- 
rection of the displacement (whether it is parallel or per- 
pendicular to the line of sight). This fact then suggests 
automatically the adoption of a cylindrical coordinate 
system for k space. Then, considering for simplicity the 
case of kikj{C a Q) we have 

hk.iCaCl) = fcf(Ca,||CM> + k\\h,x{Ca,\\C b ,x) 

+h±h(ti±tb,\\) + ki,±hx(£,±<t,±)- (B4) 

Finally, letting ki,±r l ± = k±r± cos (7) and using 
Eqs. (|A14HA20|) it is possible to recast the different terms 
appearing in the above equation as 



k }\ (Co, || Cb, || > 
k\\h,±(Ca,\\Cb,±) 



k\\ S a b, 

— kukx cos(7)V&, 



Tab 

Equation (|B2[) above then becomes 
k±dk±dkud r y 



cos(27) 



(B5) 
(B6) 

.(B7) 



(AB) = 



l_L e ^r ll+ k ± r ± co S ( 7 )] pABik ^ a ^ b) 



x exp i kS S a , 



Saa + Sbb 



2 

k\\kx cos(7) (V a - V b ) 



Dab 



cos(27) 



Tbb 



(B8) 



where for sake of clarity we have written down explicitly 
the measure of the integral and we should remind that 
since the displacement vector r makes an angle a with 
the line of sight rii = rcos(a) and r± — rsin(a). 

If we now consider the "purely angular limit" in which 
only displacements that are perpendicular to the line of 
sight are taken into account (as in (l7j ) we restrict our- 
selves to the case in which the scalar and vector terms of 
the correlators vanish (the vector term because V a =Vb, 
the scalar term because S a b — S aa — Sbb)- In this case 
Eq. (|B8|) then simplifies considerably into 



>(B9) 



x exp < kj 



T a b - ^f- cos(2 7 ) - T aa 



This expression agrees with the one obtained by Challi- 
nor and Lewis [l3| for the CMB. The point that need 
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to be stressed, however, is that it both Eq. (|B8[) and 
Eq. (|B9|) are general expressions that are valid for any 
kind of sources. Finally, the above integral can be carried 
out exactly by expanding exp(D a b cos(27)/2) in power 
series and then integrating term by term. 

Finally, it seems necessary to point out here the dif- 
ference between the two approaches and the assump- 
tions that are underlying both of them. The "non- 
perturbative" approach requires ki(Cl~Q) to be gaussian 
distributed. Once this price is paid, the apparent reward 
is to be able to take into account the full exponential, that 
is all the infinite terms appearing in its series expansion. 
On the other hand, the "perturbative" approach does 



not require such an assumption simply because higher 
order effects would contribute to higher order correlators 
and these are automatically discarded when the series ex- 
pansion of the exponential is truncated to second order. 
It seems necessary to point out, however, that the in- 
creased accuracy that can be attained adopting the first 
approach is actually hard to evaluate. If non-linearities 
are present (and this is the case when considering that the 
lensing displacements are proportional to the gradient of 
the gravitational potential, which goes non-linear at late 
epochs/low redshift) it is then questionable whether sum- 
ming all the terms appearing in the exponential would 
really lead to a consistently more accurate result. 



